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Blowup of small data solutions 
for a quasilinear wave equation 
in two space dimensions 

Q> ■ By Serge Alinhac 

o\ 

Abstract 

For the quasilinear wave equation 

dju - Au = u t uu, 

we analyze the long-time behavior of classical solutions with small (not ro- 
^ tationally invariant) data. We give a complete asymptotic expansion of the 
lifespan and describe the solution close to the blowup point. It turns out that 
this solution is a "blowup solution of cusp type," according to the terminology 
of the author [3] . 

Resume 



•3 



Pour Pequation d'onde quasi-lineaire 

d 2 u - Au = u t u u , 



nous analysons le comportement en grand temps des solutions classiques a 
donnees petites. Nous donnons un developpement asymptotique complet du 
temps de vie et decrivons la solution pres du point d'explosion. Cette solution 
est une "solution eclatee de type cusp," selon la terminologie de l'auteur [3]. 



X 



Introduction 

We consider here the quasilinear equation in R 2+1 : 
(0.1) dfu - A x u = u t u t t 

where 



xq = t, x=(xi,X2), r = \J x\ + ;c|, j;i=rcosa;, X2 = rsinu;. 
We assume that the Cauchy data are C°° and small, 

u(x, 0) = eu\{x) + e 2 u 2 ) {x) + . . . , Ut(x, 0) = eu\{x) + e 2 u\(x) + . . . 
and supported in a fixed ball of radius M. 
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Our aim is to study the existence of smooth solutions to this problem, more 
precisely the lifespan T £ of these solutions and the breakdown mechanism when 
these solutions stop being smooth. 

This problem was introduced and extensively studied by John, for this 
and more general quasilinear wave equations, in space dimensions two or three 
(see his survey paper [9] and the references therein). Then lower bounds of 
the lifespan were obtained by Klainerman ([11], [12]), Hormander ([7], [8]) and 
many other authors. Using some crude approximation by solutions of Burger's 
equation, Hormander [7] has obtained in dimensions two and three explicit 
lower bounds for the lifespan. The result for equation (0.1) in dimension two 
is 

(0.2) liminfef e 1/2 > (max (<r, a;)) -1 = f . 

Here, the "first profile" is defined as 

(0.3) rV>(<t,w) = -4= / ^=^[i?(s,u;X) - d s R(s, u, u?)]ds, 
where R(s, lo, v) denotes the Radon transform of the function v 

R(s,u>,v) = / v(x)dx. 

J XU} = S 

Hormander simply writes in his 1986 lectures on nonlinear hyperbolic equations 
[8]: 

"Even if it is hard to doubt that (0.2) always gives the precise asymptotic 
lifespan of the solutions there is no proof except that of John [10] for the 
rotationally symmetric three-dimensional case." 

In this paper, we prove Hormander's conjecture that (0.2) indeed gives 
the correct asymptotic of the lifespan. In fact, our method of proof gives 
more than that : it provides a complete description of the solution close to 
the blowup point. It turns out that the solution is a "blowup solution of cusp 
type," according to the definitions of [3]. 

Finally, to formulate more precisely Hormander's conjecture, let us intro- 
duce further useful notation and recall a previous result on upper bounds for 
the lifespan. Let u\ be the solution of the linearized problem at 0: 

d^u\ — Ami = 0, ui(x, 0) = u 1 (x), dtui(x , 0) = u\(x) . 

We have, for r — > oo, r — t > —Co, R^ being the first profile defined by (0.3), 

R<U(r-t,w) 

Ul —2 • 

Similarly, let us now define u-i by 

d^u 2 - Au 2 - d t u-yd1u-y = 0, u 2 (x,0) = u%{x), d t u 2 (x,0) = u\(x). 
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We prove in [1] that, also for r — > oo, r — t > —Co, 

for a certain smooth i?( 2 )that we call the "second profile." We assume that 
d 2 R^ has a unique positive quadratic maximum at a point (<To,wo)> an d then 
set 

fo = (^i?( 1 )(a ,a;o)r 1 , 

n = -f 2 ^i2 (2) (^o,wo). 

The result of [2] (which is also valid for general quasilinear wave equations) 
is the following. 

Asymptotic theorem (see [2]). Under the above nondegeneracy as- 
sumption on the initial data, there exists a function T e a with the following 
properties: 

i) For all N, % > T e a - e N for < e <e N , 

ii) For some C > and (Ce 2 )" 1 < t < T e a - e^, 

' < \V 2 u(.,t)\ L oo < C- 



CT£-t ~ 1 ~ T £ a - t 

The function T e a is of the form 

f e a = e- 2 {f^f{e,e 2 lne), 
where f" is a smooth function satisfying 

f e a = f + en + 0{e 2 lne). 

Thus, for numerical purposes, the asymptotic lifespan T £ a looks like the 
true lifespan T £ ; this feature would certainly make numerical experiments, 
designed to test whether or not the solution actually blows up at time T e a , 
very hard to realize. 

We prove in this work that, for equation (0.1), one has in fact T £ ~ T e a . 

I. Results and method of proof 

1. Throughout this paper, we make the following nondegeneracy assump- 
tion on the initial data. 

(ND) The function d 2 R^{a,uj) has a unique positive quadratic 
maximum at a point (ao,u>o). 

Recall that the first profile was defined in (0.3). 

For equation (0.1) with small data satisfying (ND), we have the following 
theorem. 
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Lifespan Theorem 1.1.1. The lifespan T £ of the solution u of (0.1) 
satisfies 

(1.1.1) f e = e(f £ ) 1/2 = f + efi + 0(e 2 lne). 

Moreover, for t > Tq£~ 2 (0 < to < r ) and e small, 

i) TTte solution u is of class C 1 and \u\ c i < Ce 2 ; 

ii) There is a point M £ = (m £ ,T £ ) such that, away from M £ , the solution u 
is of class C 2 with \u\ C 2 < Ce 2 there; 

hi) The solution satisfies 

(1.1.2) | V 2 u(.,t)| L oc<^^, 
(1-1-3) |5M.,t)|^>^^. 

We give here only the approximation (1.1.1) for simplicity. In fact, it is 
easily seen that the lifespan T £ and the location of the blowup point M £ can be 
computed to any order (for small enough e) by the implicit function arguments 
of [2]. In particular, T £ ~ T e a in the sense of asymptotic series. 

The inequalities (1.1.2), (1.1.3) give a rough idea of how the second order 
derivatives of the solution blow up. A much better description of the solution 
close to M £ can be obtained from the following theorem. 

Geometric Blowup Theorem 1.1.2. There exist a point M £ = 
(rh £ ,T £ ), a neighbourhood V of M £ in {(s,uj,t),s G R,uj £ S^r < f £ } and 
functions (f>,G,v G C 3 (V) with the following properties: 

i) The function (ft satisfies in V the condition 

(H) (j) s > 0,<f) s (s,LU,T) = & (s,uj,t) = M £ , 

cP ST (M £ ) < 0, V i)W (&)(M e ) = O,V^(0 S )(M £ ) » 0. 

ii) d s G = 4> s v and d s v(M £ ) / 0. // we define the map 

$(s,w,r) = (a = 4>(s,uj,t),(jj,t) 

and set $(M £ ) = (\x £ \ —T £ ,x £ \x £ \~ l ,f £ ) = M £ , condition (H) allows us to 
define near M £ a function G by 

(1.1.4) G($) = G. 

Then, close to M £ = (x £ ,T £ ), the solution u satisfies 

(1.1.5) u(x, t) = ^G(r - t, oo, et 1 ' 2 ). 
Finally, the functions <p, G and v are of class C k for e < e^. 



SOLUTIONS FOR A QUASILINEAR WAVE EQUATION 



101 



In this theorem, we see that the singularities of u come only from the 
singularities of G ; these in turn arise from the fact that the mapping <J> is 
not invertible at the point M e . More precisely, condition (H) implies that the 
singularity of <I> is a cusp singularity. Thus, describing the behavior of the 
derivatives of u near M £ is just a local geometric problem. This is the reason 
why we call this behavior of u "geometric blowup" (see [3] or [5] for details). 

2. Let us explain now the method of proof of Theorems 1.1.1 and 1.1.2. 
The idea is to construct a piece of blowup solution to (0.1) in a strip 

-C < r — t < M, r V 2 < t < T £ , < T < f 

close to the boundary of the light cone. This gives an upper bound for the 
lifespan, which turns out to be the correct one. Of course, this is not surprising, 
because the first blowup of the solution is believed to take place in such a strip, 
and not far inside the light cone. 

The proof is thus devoted to this construction, which is done in four steps, 
handled respectively in parts II, III, IV and V of the present paper. 

Step 1: Asymptotic analysis, normalization of variables and reduction to 
a local problem. We choose a number < ro < ro and use here asymptotic 
information on the behavior of u for r — t > —Co and et 1 / 2 close to ro- Thus, 
we are far away from any possible blowup at this stage, because of (0.2). 
According to [1], the solution in this domain behaves like a smooth function 
(depending smoothly also on e and e 2 lne) of the variables 

i ii 

a = r — t, lo, t = et 1 . 

Thus we set 

u(x,t) = -jj^G{a,uj,T). 

Writing equation (0.1) for G in these new variables, we are left with solving a 
local problem for G in a domain 

-C < a < M, r < r < f e , 

-1/2 

where r e = eT £ is still unknown. At this stage, we have a free boundary 
problem, the upper boundary of the domain being determined by the first 
blowup time. 

Step 2: Blowup of the problem. To solve the free boundary problem of 
Step 1, we introduce a singular (still unknown) change of variables 

$ : (s,u,t) (a = 4>(s,u,t),u,t), 4>(s,uj,t ) = s. 

The idea is to obtain G in the form 



G($) = G 
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for smooth functions <j) and G, and arrange at the same time to have (j) s vanish 
at one point M £ = (rh £ ,f £ ) of the upper boundary of the domain. Thus, we 
will have 

G s = G a (/) s , 

and the technical condition ii) of Theorem 1.1.2 gives in fact 

G a ($) = v; 

hence 

= v a /<j> a . 

We see that u, Vu will remain continuous and that V 2 u will blow up at some 
point, in accordance with the expected behavior of u. 

The nonlinear system on 4> and G corresponding to (0.1) is called the 
blowup system. 

Instead of looking for a singular solution of the normalized original equa- 
tion as in Step 1, we are now looking for a smooth solution of the blowup 
system ; however, we cannot just solve for r close to to : we have to reach out 
to attain a point where <p s = 0. 

Finally, introducing an unknown real parameter (corresponding to the 
height of the domain) , we can reduce the free boundary problem of Step 1 to 
a problem in a fixed domain. 

Step 3: Existence and tame estimates for a linear Goursat problem. Lin- 
earization of the problem obtained in Step 2 leads to a third order Goursat 
problem. In fact, it is the special structure of (0.1) which makes it possible 
to reduce the full blowup system on 4> and G to a scalar equation on 4>. The 
(unknown) point where <p s vanishes is a degeneracy point for this equation. 
Energy estimates can then be obtained using an appropriate multiplier. We 
prove in this step existence of solutions and tame estimates, which allow us to 
solve the nonlinear problem by a Nash-Moser method. 

Step 4: Back to the solution u. Having G and 0, we deduce G and thus 
obtain a piece of solution u of (0.1) with the desired properties. It remains to 
see that u = u where u is defined, and that u does not blow up anywhere else. 

II. Step 1: Asymptotic analysis, normalization of variables 
and reduction to a local problem 

1. The asymptotic analysis of (0.1) was carried out in [1]. Fix 

< T\ < T < T 2 < To- 

Introducing the variables 

1 ii 

a = r — t, to, t = et ' 
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as before, we only need here the behavior of the solution in the region 

{n <r <r 2 , r-t > -C }, 
that is, far away from any possible blowup. The result of [1] is that if we set 

u(x,t) = -^G(<t,u,t), 

the function G is bounded in C k (independently of e) for e < et (e* depends 
of course on Cq,t\ and T2). For e = 0, the function G reduces to the function, 
abusively denoted by R^\a,u),r), solution of the Cauchy problem 

(2.1.1) d T G - ^(d a G) 2 = 0, G(a,u,0) =R ( - 1) (a,u). 

According to a simple computation, the function G satisfies an equation of the 
form 

(2.1.2) -d 2 aT G + (d a G)(d 2 a G) + e 2 E(a,iu,T,G,VG,V 2 G) = 0, 

where E is a smooth function, linear in V 2 G, which we need not know explicitly. 

2. To prepare for Step 2, it is important to see that if we set w = ut 
and take the t-derivative of the equation, we obtain the conservative nonlinear 
equation 

(2.2.1) P(w) = d?w-Aw-^d 2 (w 2 ) = 0. 

Note that, with w = -rmF, 



e 2 

(2.2.2) F = dG, d = -d a + — d T . 

It 

We need the expression of P(w) in the variables a, to, r. 
Lemma II. 2. There exists the identity 



r 



(2.2.3) -P(w) = - FF aa [F aT - —F TT ] 

e z t 4r 

F 2 

- e 2 R~ 3 / 2 F^ - {F a - —F T ) 2 + e 2 hVF + e 2 h = P{F) 

It 

where h and ho are smooth functions of (to, R, r, F), and R = t 2 + e 2 a. 

We want to solve P(F) = in a (still unknown ) domain 

— Aq < a, uj G S 1 , To < r < f £ , 

with two trace conditions on {r = To} corresponding to that for u and F 
supported in {a < M} (A is big enough). 
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III. Step 2: Blowup of the problem and 
reduction to a Goursat problem on a fixed domain 

1. Formal blowup. We set, with an unknown <p, 

(3.1.1) <?($) = G, F(*) = v, $(s,u,t) = (<t>(s,u,T),u,T). 
We have then, with y = uj or r, 

(3.1.2) (d a Gm = fc'dsG, (d y G)(<f>) = d y G-fe)d s G, 

Ys 

and in particular 

(3.1.3) AG(d>) = ^AG, A = -(i + £ ^)d s + e ^d T . 

It It 

For second order derivatives of G, we find an expression of the form 



(3.1.4) (V 2 G)W = ^fA-^ + ^, 

where A and B are smooth. 



d s A B 
— 1 

<P 2 S <t>s 



Let us explain now heuristically how we establish the blowup system. 
Our aim is to substitute the expressions (3.1.2) and (3.1.4) into the equation 
(2.1.2) for G and take the coefficients of the various powers of (pj 1 to be zero. 
Of course, if we do this in a straightforward manner, we will obtain too many 
equations on G and (f>. Another possibility is to introduce an auxiliary (smooth) 
function v and force the relation 

(3.1.5) d s G = (f> s v. 

We see then from (3.1.2) that VG is smooth and V 2 G is of the form £■ + B 
(with A,B smooth); equating to zero the coefficients of 1 and of 4>~ l in the 
equation for G yields then two equations, which give, along with (3.1.5), a 
(3 x 3)-system on G,v,4>. Here, we take advantage of formula (3.1.4) and of 
the conservative character of equation (2.2.1) to get a (2 x 2)-system on v,<f>, 
as indicated in the following lemma. 

Lemma III.l. Since the functions v and F are related by (3.1.1), 
P(F)($) = -U ss v s T + -Lti + -^T 2 + T 3 , 



where 



(3.1.6) % = qv-^— <Pr-e 2 — -5-$ + 



£ 



2 



T \T 2 ^ T #3/2^' 



(3.1.7) Ti = -d s (v s %), 

(3.1.8) T 2 = Zd s v - e 2 v s N<f> + E 2 v s h 2 (ui,T,v,v LJ ,v T , . 

(3.1.9) T 3 = e 2 Nv + £ 2 h 3 (uj,T,v,v u} , v T , <f>), 
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h-2 and I13 being smooth functions and R = r + e (ft. Moreover, 

2 4 

(3.1.10) q=l + — <ft T + ^(ft 2 T , Z = 5 1 d T + e 2 5 2 d u) 



with 

(3.1.11) S 1 = --{R l/2 - e 2 v + e 2 Rl/2 ~ ^V }, 5 2 = 2R^/\ 
T 2t 



(3.1.12) N = RXl \ T f V d 2 T - R~ 3/2 dl = + 2e 2 N^d 2 TU) + N^d 2 . 



We note the three following important facts: 

(3.1.13) g^O, 

(3.1.14) <5i = -l + 0(e 2 ), 

(3.1.15) iV {1) = — + 0(e 2 ) > 0, iV (3) = --3 + 0(e 2 ) < 0. 

The fact that iV( 2 ) i s actually zero does not play a role in the subsequent 
computations, so that it is more natural to keep it. 

In order to solve the equation P(F) = 0, we now take v and <f> to solve 
the blowup system 

(3.1.16) T = 0, T 2 + S T 3 = O. 

2. Reduction to a free boundary Goursat problem. In this section, we 
are going to reduce the blowup system (3.1.16) to a scalar problem on 0, with 
boundary conditions given on characteristic boundaries of the (still unknown) 
domain. 

2.1. A local solution of the blowup system. From the implicit function 
theorem, we can write equation Tq = in the form 

(f) T = E(u,T,<t>,4> u ,v), 

with 

(3.2.1) E(oj, t, (ft, 0, 0) = 0, 

and, for e = 0, 

E = v. 

The function F being in fact known and smooth in a small strip 

Si = {t < T < T + 77 , 77 > 0}, 

we can solve, for n small enough, the Cauchy problem 

(ft T = E(u,T,4>,<t>w,F(<t>,U,T)), (ft(s,UJ,T ) = S 
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in this strip. Setting then v = F(4>,uj,t), we obtain a smooth particular 
solution (v, 4>) of (3.1.16). Note that, thanks to (3.2.1), v and 4> — s are smooth 
and flat on {s = M}. 

2.2. Straightening out of a characteristic surface. Consider the "nearly 
horizontal" surface £ through {r = tq,s = M} which is characteristic for the 
operator Zd s + e 2 cj) s N, the coefficients of Z and N being computed on (v, 4>). 
The surface £ is defined by an equation 

r = ip(s,u) + r , 

where ijj is the solution of the Cauchy problem 
(3.2.2) 

(-<5i + e 2 S 2 ^)tps + e 2 MN {1) - 2e 2 N^^ + iV (3) V£) = 0, ip(M, to) = 0. 

Equation (3.2.2) has, for small e, a smooth solution in the appropriate domain. 
This solution is 0(e 2 ) and decreasing in s. 
We now perform the change of variables 

(3.2.3) X = s, Y = w, T=(l- X (^-^))(T-To) + (T-ro-Vk( L ^), 

where x G C°°, x(t) = 1 for t < 1/2, x(t) = for t > 1. Note that this change 
reduces to T = r — To away from a neighbourhood of {r = To}. The (still 
unknown) domain 

A/, = {-^o < s < M,w G 5 1 ,r + ^ < 1" < t £ } 

is taken by this change into 

^ = {-^o < X < M,Y £ S\0 <T <T = t £ -t }. 

With a slight abuse of notation, we will again denote by (v, cj>) the local solution 
of (3.1.14) transformed by (3.2.3); this solution exists now in a small strip 
{0 < T < m} of D. 

2.3. Reduction to an equation on (f>. The equation 7o = allows us to 
express v in terms of (f> in the form 

(3.2.4) v = V{lo,tAAuAt)- 

Replacing v by V in (3.1.16), we obtain a third order equation on cp, according 
to Lemma III.l. The change of variables (3.2.3) gives 

d s = d x + T s d T = s, d U) = d Y + T ul d T , d T = T T d T , 

where 

T s = 0(e 2 ), T UJ = 0(e 2 ), T T = l + 0(e 2 ) 
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are known functions. Hence the equation on <j) becomes, in the new variables, 

(3.2.5) £(0) = (ZS)V + e 2 (S<t>)NV-£ 2 (SV)N4> + £ 2 (SV)h 2 +£ 2 (S<l>)h 3 = 0, 

where Z, N, V etc. correspond to Z, N, V etc., transformed by (3.2.3). We 
note then 

(3.2.6) Z = 6 1 dr + e?6 2 dy, 

(3.2.7) N = N^dl + 2e 2 N i - 2) d 2 , T + N<®dfy. 

Our goal is now to solve C{4>) = in D with the boundary conditions 

(3.2.8) <f>(X, Y, 0) = 4>{X, Y, 0), d T c/>(X, Y, 0) = 9 T ^(X, Y, 0), 
and (f) - X is flat on {X = M}. 

2.4. Construction of an approximate solution in the large. Note that for 
e = 0, the change (3.2.3) reduces to the translation T = t — tq, while the 
blowup system (3.1.16) is 

V = (j) T , 8 2 S V = 0. 

The initial conditions for this system are 

<KX, Y, 0) = X, d T HX, Y, 0) = -d a R^{X, Y, r ). 

Hence the value 4>o of 4> for e = is 

fo(X, Y,T)=X- TRW (X, Y, t ). 

To obtain an approximate solution valid also for large values of T, we just glue 
together the local true solution </> to 4>q: 

(x, y, t) = x (-)4>{x, y, r) + (i — y, T). 

T?l ^1 

We have then 

£(^)) = / -(°), 

where Z^ -* is smooth, flat on {X = M}, zero near {T = 0}, and zero for e = 0. 

2.5. T/te condition (H). Let us consider more closely the vanishing of </>x 
in Z). On one hand, <j>x has to vanish somewhere, otherwise the corresponding 
F and u would not have any singularities. On the other hand, as will be clear 
from the linear analysis of Chapter IV, the linearized problem corresponding 
to C(4>) =0 seems to become unstable for <px < 0. Hence we are forced to 
consider the situation where 4>x vanishes only on the upper boundary of D. 
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In a way completely analogous to what we have done in [4], we expect 4> to 
satisfy, for some point M = (m,T), the condition 

(H) X >O, cf ) x(X,Y,T) = 0^(X,Y,T)=M, 

<f> 2 XT (M)<0, V x ,y(0x)(M) = 0, V 2 XtY (<f> x )(M) » 0. 

Let us show that the approximate solution <j>^ from 2.4 satisfies this condition 
(H) at time 

(3.2.9) f = T = (max<9ii? (1) (X, Y,^)) -1 . 

Thanks to the nondegeneracy assumption (ND), the function 

d a (-d a R^\a,cj)) 

has a quadratic minimum at (<7o,u;o). On the other hand, the function 
— d a R^ (a, u>, t) is a solution of Burger's equation: at time To, its a derivative 
also has a quadratic minimum at the corresponding point, image of (oo^o) 
by the characteristic flow. In addition, To = tq — tq. Finally, dx4> > close to 
{T = 0}. 

3. Reduction to a Goursat problem on a fixed domain and condition (H). 

3.1. Reduction to a fixed domain. Recall that we want to solve the equa- 
tion C(4>) = in a domain such that 4> satisfies the condition (H) for a point lo- 
cated on the upper boundary. The approximate solution , starting point of 
some approximation process, satisfies this condition for a domain of height To, 
according to 2.4, 2.5. Unfortunately, in the successive approximation process, 
further modifications of (j)^ will yield functions not satisfying (H) anymore. 
We are thus forced, at each step of the process, to adjust the domain to have 
the new 4> satisfy condition (H). 

To achieve this, we introduce a real parameter A close to zero, and perform 
the change of variables 

(3.3.1) X = x, Y = y, T = T(t, A) = T (t + Xt(l — Xi (*)))> 

where xi is 1 near and near 1, and To is defined as in (3.2.9). Of course, 
one should not confuse these variables with the original variables ! We will 
from now on work on a fixed domain 

D = {-A < x < M, yeS 1 , < t < 1}. 

We denote the transformed equation by 

(3.3.2) C(X,(f>) = 0, 

the transformed approximate solution for A = A*- ) = by 

(3.3.3) <A (0) (x, y, t) = (x, y, T(t, 0)) = ^ (x, y, T t), 
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and set 

(3.3.4) C(X^,^) = f^ = P\x,y,T t). 

We note that <f>^ satisfies (H) in D$ for a certain point Mq = (mo, 1). 

3.2. Structure of the linearized operator. The linearized operator of £ at 
the point (A, (j>) is denoted by 

(3.3.5) £'(X, 0)(A, $ = d\£(X, <j>)\ + d^£(X, <j>)j>. 

Because £(X,(f>) comes from £(<p) by (3.2.1), we have the following lemma. 
Lemma III.3.1. If £(X, </>) = f, then 

(3.3.6) d x £(X, <j>) + d^X, 4>)(dt4>^) = d t f^. 

For the time being, it is not necessary to make an explicit computation of 
dffrC. Note only that, if we have at some stage £(X,(p) = f (for a small /), to 
solve 

£'{X,<p)(\4>)=f 
approximately it is enough to solve 

d 4> £(x,4>)* = f 

and then to take (A, (ft), verifying 
In fact, we get with this choice 

? d x T 



£'{\4>){XA) = f-\d t f 



The additionnal term contains a product of small A by small dtf, which is 
negligible as a quadratic error. Having determined we see that we still have 
an additional degree of freedom to choose 4>: we will take advantage of this to 
arrange for (f> + <fi to satisfy (H) . 

3.3. The fundamental lemma. We follow here exactly the same idea as in 

[4]- 

Lemma III. 3. 2. Assume that (fr—cf)^ and ip are small enough in C 4 (Do). 
Then 

i) If 4> satisfies, for a certain m, 

4>x(m,l) = 0, V Xj y(</) x )(rh,l) = 0, 
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it also satisfies (H). 

ii) There exist a function A((p, ip) and a point fh((f>, ip) such that 

A(^(°),0) =0, m(<^°\0) =rh 

and the function <p + ip + df(pA((p, ip) satisfies (H) in Dq for the point 
M = (m(^),l). 

iii) 7/0 already satisfies (H) /or a pom£ M = (rh, 1) c/ose to M , i/ien 

A(0,O) = 0, m = m(0,O). 

Proof. Point i) is clear from the Taylor expansion. Let now 

G:(<p,iP,rh,\)^(d x <!>(rh,l), ^$(m,l), ^$(m, 1)) = (G 1: G 2 , G 3 ) 

with $ = <p + ip + \d t <p. The function G is of class G 1 from G 3 x G 3 x i? 3 to 
i? 3 . By construction of (p^\ 

G(^°\0,m ,0) = 0. 

On the other hand, 

9 A Gi(^ ),0,mo,0) = ^°)(mo,l), 

^C?i(^ ),0,trio,0) = (0,0), 

d m (G 2 , G 3 )(^°\ 0, m , 0) = V 2 (#)(m , 1) » 0. 

Hence the implicit function theorem yields A = A(<p, ip) and rh = rh((p, ip) with 
the desired properties. Thanks to i), <p + V* + dt(pA(<p,ip) satisfies (H). 
Finally, under the assumptions of iii), G(cp, 0, rh, 0) = 0; hence 

A(0,0) = 0, rh(4>,0) = rh. 

3.4. Back to the linearized operator. We go back to Section 3.2 and explain 
now how we can solve the linearized operator and get cp + <p to satisfy (H). 
Assume that <p already satisfies (H) for rh close to rho and \<p — <P^\c 4 (D ) 
small. We will have 

(3.3.7) ^ + < j ) = (j) + q f + xd t( p^. 
We now take 

(3.3.8) A = (1 + A)A(</>, 4>) = (1 + A)(A(& ) - A(0, 0)). 

Because / is small, and A are also small: the right-hand side (p of (3.3.7) is 
then close to </>(°) and satisfies at rh = rh((f), 4>) 

4> x {m,l) = 0, V(4)(m,l) = 0. 

According to point i) of the lemma, (p satisfies (H) . 
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4. An iteration scheme for the problem. To solve the problem C(X, 4>) = 
in Do, we will use a Nash-Moser scheme. We refer to [6] for notation and 
details, and specify here only the nonstandard points. 

4.1. Spaces and smoothing operators. We will work with the usual Sobolev 
spaces H s {Dq). In the process of solving, we note that our starting function 
</>(°) satisfies the good boundary conditions, so that all the modifications (f> 
we will have to consider will be "flat" on {t = 0} and {x = 0}. Hence the 
smoothing operators used have to respect this "flatness". To achieve this, we 
take a smooth function ip supported in {t > 0, x < 0} whose Fourier transform 
vanishes at the origin of order k. Setting 

S e = tp e -i*, i> £ = e" 3 V'(xe" 1 ,ye~ 1 ,te~ 1 ), 

we see that the operators Sg satisfy the usual properties: 
i) \Sgu\ s < C\u\ s >, s < s', 
h) \S e u\ s < C6 s ~ s '\u\ s/ , s > s', 

iii) \u — Sgu\ s < C9 s ~ s '\u\ s i, s < s', 

iv) \4qSqu\ s < C6 S ~ S _1 |ti| s ', but only for < s,s' < s^, s^ going to infinity 
when k goes to infinity. Here, |.| s denotes the i/ s -norm in Dq. 

4.2. Smoothing of <f>. In a Nash-Moser procedure, instead of solving 
at a given step, we solve 

d <t> £(\,4>)v = f 

for an appropriate smoothing <fi of <fi. Wishing to have <p satisfy condition (H), 
we use the following twin of the fundamental Lemma III. 3. 2. 

Lemma III. 3. 2'. Assume \<p — 0^|c 4 (-D o ) small enough. Then there exist 
functions A.(<p) and rh{<f>) such that 

A(4> {0) ) = 0, m(^°)) = m 

and 4> + xA(4>) satisfies (H) for the point M = (rh, 1). Moreover, if cp satisfies 
(H) for rh close enough to rho, 

A(4>) = 0, rh(4>) = rh. 

4.3. Approximation scheme. Assuming that we can solve the equation 
d^C^i = / in flat functions, with a tame estimate (see Propositions IV. 3. 2 and 
IV. 4 for a precise statement), we set 

e n = (9 r f+n)\ S n = S dn , 

0(n+l) = ^(n) + A ^(n) ) A (n+1) = A (n) + AA (n)_ 
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Here, A means "modification" and has nothing to do with the Laplacian! The 
parameters 9q and 77" 1 are chosen big enough. For the special smoothing of 
(p^ discussed in 4.2, we set 

S n <j>W = S n <f>W + *A(S n (n) ) = + x(A(S n ^) - A(0("))). 

Knowing \^ n \^- n \ we solve in flat functions 

for 7„ to be determined. Then we take 

AAW = (l + A( ri ))[A(0("),^ n )-A( ( /)( n ),O)], 
A^W = * n + AA^a^W^(AW,t). 

We now determine the 7n . First we define the three errors of the solving 
process: 

i) The Taylor error is 

4 = £(A( fc+1 \ 0( fe+1 )) - £(AW, - £'(A«, 0«)(AA« A</> (fc) ). 

ii) The substitution error is 

4' = [fy£(AW ^< fc >) - ^(A( fc ),S fc ^ fc )]* fc . 

iii) The result error is 

e' fc " = AA( fc )||(t,A( fc ))a t (£(AW,#))). 
Then we see that 

£( A (n+D )(/) (n+l)) _ £( A H, = e n + 7n , 

where e n = + e" + e"' is the total error. Finally, we denote by 

E n = So<fc<ri-iefc 

the accumulated error, and 

£(A^ +1 \ 4> {n+l) ) = / {0) + S < fc < n7fc + E n+1 . 
It is natural at this stage to determine the 7n by 

So<fc<n7fe + S n E n = —S n f(°\ 

which leads to 

£( A ("+ 1 ), = /(°) - 5n /(0) + e „ + (£„ - S n £ n ). 
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IV. Existence and tame estimates for the linearized problem 

1. Structure of the linearized operator. In part III, we showed how solving 
£(A, <p) = can be reduced to solving the linearized equation d^C. We display 
now the structure of this operator. 

Proposition IV. 1. The linearized operator has the form 
(4.1.1) cofy£(A,0) = ZSZ + e 2 {S(/))NZ + e 2 Bd 2 + b SZ + e 2 £. 

Here 

i) Z is a field of the form 

Z = d t + e 2 z d y ,z = zo(x,y,t,X,^,<j>y,^ t ), 

ii) S is the field (independent of </>) 

S = d x + e 2 s dt, s = s (x,y,t, A), 
hi) N is the second order operator 

N = N X Z 2 + 2e 2 N 2 Zd y + N 3 d 2 , N t = N^x, y, t, A, </>, fa, <j> t ), 

*■ = + oi£2) - * 3 = irJnt,X )f + 

iv) B = -(ZS<P)N 3 , 

v) c = c (x,y,t,X,4>,(j) y ,(j) t ) = (qd t T)~ 2 + 0(e 2 ), 
b = b (x,y,t,\,V(f),V 2 (i)), 

vi) £ is a second order operator which can be written as a linear combination 
of 

id, S,Z,d y ,SZ,Z 2 ,Zd y ,d 2 

with coefficients depending on the derivatives of 4> up to order 3. Moreover, 
£ does not contain d 2 for e = 0. 

Proof, a. The linearized operator d^C is obtained as follows: first, we 
linearize the 4> equation resulting from substituting (3.2.4) into (3.1.16). Then 
we perform the changes of variables (3.2.3) and (3.3.1). 

b. With the notation of Lemma III. 1, let us compute Tq. We find 

T = qv + Zip + e 2 7i0, 71 = 71 (w, r, 4>, 4> W ,4> T ). 

Hence V = —^(Z<p + e 2 ji4>). On the other hand, linearizing T 2 + <f> s Tz = 
gives 

ZV S + ZV S + e 2 <p s NV + £ 2 <f) s NV + e 2 (f> s NV - e 2 V s N<j) 

- e 2 V s N<j) - e 2 V s N<j) + e 2 V s h 2 + e 2 V s h 2 + e 2 <j) s h 3 + e 2 <j) s h 3 . 
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We see that Z and N yield only first order derivatives of <fi, multiplied by e 2 . 
The same is true of hi and /13, except for the terms V w and V T coming from 
the corresponding terms v u and v T in h 2 and It follows that the linearized 
equation on (j) has the form 

{Zd s + e 2 <f> s N - e 2 {N<p)d s + e 2 h 2 d s + e 2 h 4 d UJ 

+ e 2 h 5 d T }V - e 2 V s N<p + e 2 x V0, 

the last term denoting just first order derivatives of </>. 

c. The composition of the two changes of variables, denoted by bars, 
operates the following transformation of operators (to avoid introducing un- 
necessary notation, we denote by * known functions): 

d 8 = d x + e 2 s (x, y, t, X)d t = S, 

B UJ = dy+e 2 * {x,y,t,X)d t , 

3 T = (l + e 2 *(x,y,t,\))(d t T)- 1 d t , 

Z = (-1 + s 2 * (x, y, t, A, 4>, 4> y , 4> t ))Z, 

N = N {1) dl + 2e 2 N i2) d 2 yt + N (3) d 2 + lower order terms, 

= 4( dt T) 2 (rl + T(t,X)) + ^ * (X ' y ' A ' +> *vM> 

^ (3) = ~ ( TQ + T(t A)) 3 +g2 * 

Finally, we replace the operators by their transforms into the linearized equa- 
tion and set N = —(dtT)N. To obtain the value of B, we observe, keeping 
only the e 2 terms in the coefficient of d 2 , that 

B = q(d t T) 2 (SV)N^ + 0{e 2 ), 
from which iv) follows. □ 

2. Energy inequality for the linearized operator. In the following, in a 
(desperate) attempt to simplify the notation, we will write abusively Z for Z, 
S for S, N for N, and replace e 2 everywhere by e. 

We assume a given smooth function (j) in Dq , close to <fi(°) , satisfying (H) 
for a point M = (m, 1). We then set 

P = ZSZ + e(S<j>)NZ + eBd 2 + et + b SZ, 

P = ZSZ + e(S<f>)NZ. 
Recall that we have arranged for {t = 0} to be characteristic, that is 
(4.2.1) (t = 0) s + (S<f>)Ni = 0. 

The connection between P and P is explained in the following straightforward 
lemma. 
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Lemma IV. 2. With the above notation, 

ZPu = PZu+£lZu+eZlu+Z(boSZu)+e[ZB-2£B(dyZ )]dyU-e 2 B(dyZ )dyU. 

Here, £ is a second order operator which can be written as a linear combination 
of Z 2 , Zd y , d 2 , d y , the coefficient of d 2 being a multiple of(S<fi). 

The point of this lemma is that P does not contain the delicate term 
sBd 2 , as P does. 

2.1. Energy inequality for P. The following computation is very close to 
that of [4, Section III. 2]. Unfortunately, N is different here, causing some new 
problems. Thus we think it better to give the whole computation again, which 
is, after all, the heart of the proof. We set 

A = S(p, 5 = 1 — t, g = exp h(x — t), 
p = Si exp-(x-t), | - 1 o = I-Il2(d )- 

Proposition IV.2.1. Fix \i > 1. Then there exist C > 0, e > 0, r/ > 
and ho such that, for all 4> satisfying (H), (4.2.1) and \(f> — (f>^\c A {D ) ^ Vo, 
< e < eq, h > ho and smooth u with 

u(x, y, 0) = d t u(x, y, 0) = 0, u(M, y, t) = 0, 

the following inequality exists: 
(4.2.2) 

h\pSZu\l + h\pZ 2 u\l + £h\pd y Zu\l + £ 2 J 5^ l g{S<j)){l + 5h)\d 2 yu\ 2 < C\pPu\l 

Proof a. With a still unknown multiplier 

Mu = aSZu + £cd 2 u + dZ 2 u, 
we find by the usual integrations by parts 

J PuMudxdydt 

= j Ko{SZu) 2 + j£ 2 K 1 (d 2 u) 2 + j £K 2 {d y Zu) 2 + j K 3 (Z 2 u) 2 

+ £ J (dyZu)(SZu)[dy(c - aAN 3 ) - £Z(aAN 2 ) - 2£ 2 aAN 2 dyZ ] 

+ £ j \d 2 u)(SZu)[-Zc + £cd y z ] 

(Z 2 u)(dyZu)[-dy(dAN 3 ) + £S(aAN 2 ) + edyicAN^ 
+ ea^iViSzo — N 3 d y so — £ 2 zod y so) + £cd y so — 2£ 2 dAN 2 d y zo] 



+ e 3 
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+ e 2 J {Z 2 u){d 2 y u)[-Z{cAN 1 ) + ecANrfyZo] 
j {d 2 y u){d y Zu)[-Z{cAN 2 ) + ecAN 2 d y z ] 
: J (Z 2 u)(SZu)[-Z(aAN 1 ) - ed y {aAN 2 ) - eaANi 

+ e 2 j (d y u)(d 2 z )[cSZu + cAN 1 Z 2 u + cAN 3 d 2 u + ecAN 2 d y Zu] 

+ h-h + J 1 . 
The coefficients K{ of the quadratic terms are 
2K§ = — Za — eadyZo, 
2Ki = - Z(cAN 3 ) + 3ecAN 3 d y z , 

2K 2 = S(aAN 3 - c) + Z(dAN 3 ) - eZ{cANi) + e 2 d y (cAN 2 ) - e 2 dAN 3 
+ E 2 aAN 3 (Zso + ezodySo) — ce 2 (Zsq + sz^dySo — eANid y zo) 
+ 2e 3 aAN 2 (Sz + ez Zs ), 
2K 3 = -Sd + eSiaANt) - eZ(dANx) 

- ed(d t s + eANidyZo) - 2e 2 d y {dAN 2 ) 
+ E 2 aA(N 1 d t s - 2NiZs G - 2eN 2 d y s ). 
The terms Iq, h, J\ are boundary terms. We have 

2/i= / [a(SZu) 2 + edAN l (Z 2 u) 2 
J{t=i} 

- e(dAN 3 - EcANi)(d y Zu) 2 + e 2 cAN 3 (d 2 u) 2 

+ 2e(SZu)(aAN 1 (Z 2 u) + eaAN 2 d y Zu + cd 2 u) 

+ 2e 2 cANi(Z 2 u)(d 2 u) + 2e 3 cAN 2 (d y Zu)(d 2 u)]. 

The Iq term on {t = 0} vanishes, thanks to (4.2.1) and to the fact that the 
first two traces of u are zero. Finaliy, 

2Ji = / {-d + eaANi){Z 2 u) 2 - e(c - aAiVg)^^) 2 



L 



x=-A } 

+ 2e 2 aAN 2 (d y Zu)(Z 2 u). 

b. We choose now 

a = A~ 1 5^g, c = c'5^g, d = -d'5^g 

where h,c',d' are positive constants to be chosen later. We analyze first the 
boundary terms. 

We find that 2I\ is the integral of 

a(SZu + eAN 1 Z 2 u + e 2 AN 2 d y Zu + c'Ad 2 u) 2 
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- e(Z 2 u) 2 5^gAN 1 (d' + eiVi) + e{d y Zu) 2 8^gA{dl ' N 3 + ec'TVi - e 3 iV 2 2 ) 

+ e 2 (d 2 u) 2 5»gc'A{N 3 - c') - 2e 3 A5^gN 1 N 2 (Z 2 u){d y Zu). 

From Proposition IV. 1, it is clear that I\ > for d < N 3 and e small enough. 
The term 2J\ is the integral of 

(5"g)[{d' + eiVi)(Z 2 n) 2 + £(iV 3 - d){d y Zu) 2 + 2e 2 iV 2 (5j / Zn)(Z 2 u)]. 

Clearly, Ji > for d < N 3 and £ small enough. 

c. We analyze now the signs of the quadratic terms. 
We find 

2K = —^-\bZA + fiA + 5A(h - ed y z )}. 
In a small neighbourhood of M, we have A = <p x > — 8V, with 

^<0, y(M) = ^(M)<0 
according to (H). It follows that, with // = ^±1 and a possibly smaller w, 
5ZA + //A > (5(^ t + ez <t> xy - n'V) > 0; 

hence 

#o > C^-^(l + 5/i). 
Outside w, we have for h big enough 

<5(ZA + A(h - ed y z )) > 1/25 Ah; 

so finally 

^o>^o— + ^o>0. 

Now 

2#i = c'd^gANsifi + 36(5^20 ) + c'S^g[(hN 3 - ZN 3 )A - [ZA)N 3 ]. 
Close to M, < 0, hence 

(/iiV 3 - ZiV 3 )^ - (ZA)N 3 > ChA. 
This also holds away from M for big h. Finally, 

K x > k lC !5^- l gA{l + 5h). 
The first two terms in are 

5»g{SN 3 ) + g(N 3 - c')(h5»(l - es ) - e^" 1 ) 



+ nS^d'gANs + d'S^g[{hN 3 - ZN 3 )A - (ZA)N;_ 



3 ■ 



For d < N 3 , they are bigger than 

2k 2 g{5' 1 h + 5^- 1 A). 
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All the other terms in K 2 are bounded by 

CeS^gA + CehS^g. 

Hence 

K 2 > k 2 5 fl - 1 g(A + 5h). 

Finally, 

2K 3 = g{d' + eiV 1 )[M^(l - es ) - e^qS^ 1 } 
+ eS^giSNr) - eN 1( f 'AS"' 1 g(ji + hS) 
+ 0(e5»g) > k 3 g(h5» + ed^A). 

d. Consider now the product term in (d 2 u)(SZu). To ensure positivity, 
it is enough to check the positivity of the quadratic form 

-9Za(SZuf - e 2 eZ{cANz){d 2 y u) 2 - 2e(Zc)(d 2 u)(SZu) 

for some 9 < 1. The discriminant A satisfies 

-e~ 2 A = e 2 (Za)(Z{cAN 3 )) - {Zcf. 

After rearrangement of terms we get 

_4<5 2 A 2 e- 2 A = (6 2 N 3 - c')A 2 (fi + hS) 2 - 6 2 N 3 5 2 (ZA) 2 

- 6 2 A 2 {ZN 3 )5(fi + hS) - 9 2 5 2 A(ZA)(ZN 3 ). 

Away from M, A is uniformly positive and the right-hand side is positive for 
big h if d < 9 2 N 3 . Close to M, 

A > -SV, V(M) = ZA- ez d y A, 

and the right-hand side is bigger than 

{(9 2 N 3 - cV + 0(5)}5 2 V 2 - 9 2 N 3 5 2 (ZA) 2 . 

For d small enough and 9 close to 1, this is bigger than C5 2 . 

e. It is easy (but lengthy) to check that all other product terms can be 
absorbed in the quadratic terms for small e. 

f. Finally, using Cauchy-Schwarz inequality, we obtain (4.2.2). □ 

2.2. Energy inequality for P. We deduce from (4.2.2) and Lemma IV. 2 
an energy inequality for P. 

PROPOSITION IV. 2. 2. Fix n > 1. Then there exist C > 0, e > 0, rj > 
and ho such that, for all <p satisfying (H), (4.2.1) and \<p — (j>^\c 4 (D ) — Vo, 
< e < £q, h > ho and smooth u with 

u(x,y,Q)=d t u(x,y,0)=0, u(M,y,t)=0, Pu(x,y,0)=0, 
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there exists the energy inequality 

(4.2.4) h\pSZ 2 u\l + h\pZ 3 u\l + eh\pZ 2 d y u\l + e 2 \ P Zd 2 u\l 

+ e 2 h J 5^g{S(f>)\d 2 yZu\ 2 < C\pZPu\\. 

Proof, a. First, let us compute in Pu the coefficient of d 2 u on {t = 0}. 
We get 

[d x + e 2 (2s z + 2AN 2 + 3Az )d y + e 2 z dyS ](d?u). 

Hence the assumptions on u imply d 2 u(x, y, 0) = 0. 

b. We can then apply inequality (4.2.2) to Zu, getting 

h\pSZ 2 u\ 2 + h\pZ 3 u\l + eh\pd y Z 2 u\l 

+ e 2 j S^giS^il + 5h)\d 2 yZu\ 2 < C\pPZu\\. 

We now have to bound the terms of PZu — ZPu (given in Lemma IV. 2) by 
the left-hand side. Note first that, by standard lemmas (see for instance [2, 
Lemma 2.1]), the terms of the left-hand side also give control of 

eh\pZ 2 d y u\l + eh\pZd y Zu\l + h\pZSZu\l. 

Thus the terms from e(£ + £)Zu are easily dominated for small e and big h. 

Since 5 < Ccf) x , we have control of e 2 \pd 2 Zu\Q, implying control of both 
e 2 \pZdyU\ 2 j and e 2 h\pd 2 u\ 2 ] . This makes it possible to absorb the terms e[Z,£] 
and also the remaining terms in d 2 u, d y u, SZu and ZSZu for small e and 
big h. □ 

3. Higher order inequalities. 

3.1. The spaces H s . In Section 2, we used the fields Z, S, d y systematically 
instead of the standard set dt,d x ,d y . The reason for this is that if we develop 
the expression of P, the energy inequality becomes much less transparent (to 
say the least) ; in fact, one should observe that (4.2.2) by itself does not give 
separate control of d 2 y u, for instance. Thus it is also appropriate to commute 
P with Z, S, d y and their products. We define first T l as any product of I fields 
among Z, S, d y , and for integer s, 

H s = {n G L 2 (D ), l<s^T l ue L 2 (D )}. 



We denote the natural norm by 
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Proposition IV. 3.1. The spaces H s and H s are the same. Moreover, 
for <f> bounded in C 1 , 

i) Hl<c(H2 + Hi=o(i + H2 +1 )), 

ii) H2<c(Hl + Hi=c(i + H5 + i)). 

Proof, a. The equality of the spaces is obvious for smooth (j). 

b. We will repeatedly use the following classical lemma (see for instance 

[6])- 

Lemma. For \a\\ + . . . + \a p \ = s, 

\d ai ui . . . d ap u p \ < CEi<j<p|iii|L°° • • • \uj\ s . . . \u p \ L ™. 

Denoting by d usual derivatives, we have (skipping everywhere irrelevant 
numerical coefficients) 

T l u = Zd qi a...d qk ad r u, 

where a = eso or ezq, and 

r > 1, k < I, \qi\ + . . . + \qk\ + \r\ = I. 

Hence 

\T l u\ < C(\u\i + \u\ L oo\a\i), 

which gives i). 

c. Conversely, we have 

d t = Z-ez d y , d x = S - es Z + e 2 s z d y . 
For p + q + m = Z, we can write 

dfdZd™ = ST 91 (a) . . . T qk (a)T r u, 

with 

\qi\ + . . . + \q k \ + \r\ = p + q + m, k < p + q + m 

and a = ezo, esq or e 2 soZo- Using the identity of b, except for the terms T l u, 
we get 

dldld™ = So . . . aT l u + Ed Pl a . . . d p ^ad r u, 

with \pi\ + • • • + \pj\ + \r\ = I, r < I — 1 in the last sum. Taking L? norms and 
using the interpolation lemma, one obtains 

|fiffiga^«| < C|T^|o + C|«|,_i + C|u| L oc|a|,. 
Hence, by induction on s, we have ii). □ 
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3.2. Structure of the commutators. The two following straightforward 
lemmas describe the structure of the commutators of P with a product K = T l . 

Lemma IV. 3.1. Let us denote here by c various coefficients which are 
smooth functions of 

x,y,t,\,4>,V<t>, V 2 4>, 

and by Q3 the principal part of the third order operator Q . Then 

i) [Z, P] 3 = e(cZ 3 + cd y Z 2 + c&Z), 

ii) [S, P] 3 = e(cZ 3 + cdyZ 2 + cd\z + cZSZ + cd y SZ + c(S(t>)d 3 ), 

iii) [d y ,P] 3 = e(cZ 3 + cd y Z 2 + cd 2 Z + c(S<f>)d 3 ). 

Lemma IV.3.2. When K = T l , 
Z[K, P]u = ZX[T, P\zT l - l u + ZY>d qi a . . . d qk ad q bT r u 

+ £d 9l a . . . d qk ad q bT l u + Zd qi a . . . d qk ad q bd r u 
= ZSi + ZS 2 + S 3 + S 4 . 

Here, a = eso or ezq, b is a coefficient of P, and the conditions on the deriva- 
tives are the following: 

i) In £ 2 , i < r < I + 1, \q\ + \qi\ + ■■■ + Wk\ + \r\ < I + 3. 

ii) In £ 3 , |g| + + . • • + |%| < 4. 

iii) In £4, r <l- 1, |g| + + . . . + |qr fc | + |r| < / + 4. 

Moreover, all the terms in Z[K,P]u contain a factor e, except for 
Z[K, boSZ]u which is of the special form 

Z[K, b SZ]u = Z 2 £ . . . T l u + ZSYi . . .T l u + Z£ . . .T l u 
+ £ . . . T l u + £ r </_i . . . d r u. 

3.3. Higher order tame estimates. 

Proposition IV.3.1. Let (j) satisfy the assumptions of Proposition 
IV. 2.1 and moreover \4> — <t>^\c 6 (D ) — Vo- Then there exists Eo such that for 
< e < eo and all s, there exists C s such that for all smooth u flat on {t = 0} 
and {x = M}, 

(4.3.1) \u\ a < C\d^Cu\ s+1 + C|fy£u| 3 (l + \<f>U+7). 

Proof, a. Recall that c^d^C = P. Thus, it is enough to prove (4.3.1) for 
P. Let Pu = f,K = T l ,l < s. Then 

KZf = [K, Z]f + Z[K, P]u + ZP(Ku). 

Applying the energy inequality (4.2.4) to Ku and summing over I < s, we 
obtain control of a certain number of terms which we denote by E s and do not 
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repeat here. We have to show that the terms of Z[K,P]u, whose structure is 
discussed in Lemmas IV. 3.1 and IV. 3. 2, can be bounded in the weighted norm 
\p.\o by an arbitrarily small fraction of E s (in short, "absorbed"). 

b. Examine first the terms of ZH\. Taking into account the explicit 
expressions of the commutators [T, P]% given in the lemma, we see that there 
are six different types of terms to control, which we display in three groups: 

(1) eZ(cZ 3 )T'- 1 n, eZ{cd y Z 2 )T l - x u, eZ{cZSZ)T l ^u, 

(2) eZ(cd 2 y Z)T l - x u, eZ{cd y SZ)T l - x u, 

(3) eZ{c{S4>)dl)T l - x u. 

The terms of group (1) are clearly absorbed for big h. 
For the group (2), we write 

dlZT l - l u = dyZ(d y T l - l u) + *dy{d y T l - l u) + *(d y T l - 1 u), 

which shows absorption for big h. We proceed analogously by splitting the 
second term in terms easily absorbed 

d y SZT l - l u = SZ{d y T l - l u) + *S{d y T l ~ 1 u) 

+ *Z 2 T [ - l u + *d y ZT l - x u + *d y T l - x u. 

Finally, 

z( c {s<p)dl){d y T l ~ l u) = (Zc){s<p)d 2 y {d y T l ~ l u) + cZ{S4>)d 2 y {d y T l ' l u) 

+ c{S<P)d 2 y Z(d y T l - l u) + c(S<f>)[Z, d 2 y\{d y T l - l u). 

Since [Z, d 2 ] = *d 2 + *d y , all the terms can be absorbed for big h. 

c. We analyze now ZS2 and S3. 

A product T s+1 containing at least one factor S can be written 
T s+1 u = ST s u + Y,daT s u + £<9 91 a . . . d qk ad r u, 
r < s - 1, |r| + \qi\ + . . . + \qt\ = s + 1. 

We can proceed similarly if T s+1 contains at least a factor Z. In all cases, 

eh\pT s+l u\l < CE S + Ce/i|pS r < s _i . . . d r u\l. 

A similar analysis gives the same estimate for terms like eh\pZT s+l u\^. Thus, 
all terms from ZE2 or S3 which have e as a factor are absorbed for big h, 
modulo an additional term e/i|pS r < s _i . . . d r u\Q on the right-hand side of the 
inequality. The same is true for the terms from Z[K, b^SZju, thanks to their 
special structure. 

d. We now fix h and use the inequalities 

\pv\o < C\v\ , E s > Ch\u\ s . 
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Applying the interpolation lemma for the index s — 1 ( "taking out" five deriva- 
tives from the coefficients), we can bound the additional terms and the terms 
of £ 4 : 

\d qi a...d qk ad q bd r u\ < C\u\ s -i + C|u| L oo(l + \<f>\ s+7 ). 
By induction on s, using Proposition IV. 3.1, we finally obtain 

\u\ 2 s < C\f\ 2 s+1 + C|/H(l + |<^ +6 ) + C\u\l{\ + \<f>\ 2 s+7 ). 
Since it is easy to obtain a low norm estimate 

M 2 <c|/| 3 , 

we get (4.3.1). □ 

4. Existence of flat solutions. 

Proposition IV. 4. Let 4> satisfy the assumptions of Proposition IV. 2. 2. 
Then there exists, for all f £ C co {Dq), flat on {t = 0} and {x = M}, a unique 
smooth solution of 

d^Cu = f, u(x, y, 0) = d t u(x, y, 0) = 0, u(M, y, t) = 0, 

satisfying the estimates (4.3.1). 

Proof, a. Define, as in [4], the smoothed operators P a by replacing d y by 
Y<x = x( a D y )dy in P (% being compactly supported and one near zero). For 
fixed a > 0, we can solve 

P a u = f, u(x,y,0)=d t u(x,y,0) = 0, u(M,y,t) = 

in smooth functions. In fact, setting 

S = S + eiSflmdt = d x + e(s + (SflN^dt = d x + es d t , 

we can expand the terms of P a and write 

P a u = dtSu + d t SAiu + d^A 2 u + d t A 3 u + SA 4 u + A 5 u. 

Here, the Ai are zero order operators in y depending smoothly on (x,t). To 
solve in Dq, we proceed as usual by writing the equation in integral form. 
More precisely, let Mq = (xo,yo,to) be some point in Dq and denote by 

X(M ,x),Y(M ,x),T(M ,x) 

the parametrization by x of the segment of the integral curve of S from Mq. 
We set 

(I § v)(M )= / v(X,Y,T)dx, 

J M 
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so IgV = for x = M and Slgv = v. The key point is that this segment of 
integral curve (for xq < x < M) is well defined and contained in Dq because 
of the following two facts: 

i) On {t = 0}, S = 8 X , 

ii) Since so = close to {t = 1} and 4> satisfies (H), the integral curves issued 
from points Mq where t = 1 enter Dq for x > xq. 

We similarly set 

(Iv)(x,y,t) = / v(x,y,s)ds 
Jo 

and write the equation P a u = f, 

u + IB lU + IgB 2 u + IgIB 3 u + I 2 B 4 u + I § I 2 B 5 u = I § I 2 f 

for appropriate zero order operators Bi in y that we now determine. Applying 
d 2 S to the left, we get 

d?Su + d?SIB lU + d 2 t B 2 u + d t B 3 u + d?SI 2 B 4 u + B 5 u = f. 

But 

d 2 SI = ed t (d t s )+d t S, 
dfSI 2 = 2sd?s I + 2ed t s + S. 

Hence we can take 

B x = A U B 2 = A 2 ,B 4 = A 4 ,B 3 +ed t S B 1 = A 3 ,B b + 2e(d t s G + d 2 t s G I)B 4 = A 5 . 

A standard fixed-point argument yields the smooth flat solution u a . 

b. One can prove, exactly as we have done for P, higher order energy 
estimates for P a with constants independent of a (see [4]). We can find some 
subsequence of u a converging, say, weakly in L 2 , to some solution u of Pu = /. 
Since \u a \ s is bounded for all s, u is smooth and has zero traces. It is impor- 
tant here that, according to Proposition IV. 2. 2, we do not have to decrease e 
with s. □ 

V. Going back to the solution 

1. In parts III and IV, we obtained a solution (v, 4>) of (3.1.16) in a domain 

D = {-A <s<M, lo e S 1 ,t <t < f £ }. 

If < e < £fc, this solution is, say, of class C fc .We define now a function G 
supported in {s < M} by 

ClG = (j) s V, G(s,lo,t ) = G(s,(jJ,T ). 

Note that 
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Thanks to condition (H), we see that (3.1.1) yields continuous functions 
G and F defined in 

$(£>) = {(/>(- A ,u,t) <a<M, weS 1 , r < r < f e }. 

Moreover, F satisfies there equation (2.2.3) P(F) = 0. Of course, F and G 
are of class C k for r < f £ , and VF becomes infinite at the point $(M £ ) of 
the upper boundary of ^(-D). Returning to the original variables and defining 
u = G, we obtain u and w in 

Six = {0(-A o ,cj,et 1/2 ) <r-t<M, r < rf 1/2 < f e }. 

Moreover, -u; satisfies P(u;) = there, with Vw becoming infinite at a point 
M £ . Finally, for C to be determined, define 

J) = {-,4 + C(e 2 t - r 2 ) < r - t < M, r < et 1/2 < f £ } C fii. 

We take C big enough to have 

i) t — Aq + C(e 2 i — Tq) — (t + 0(— ^4o, e ^ 2 )) an increasing function, 

ii) an influence domain for the operator d 2 — A — wd 2 . 

In this process, we can arrange for to contain M e , decreasing e if nec- 
essary (because (j) is arbitrarily close to the solution for e = 0, which we can 
arrange). The function u satisfies the equations 

d t u = w, (d 2 - A)u - (d t u){d?u) = 

in £1 and has the same traces as u on {et 1 ^ 2 = To}. By uniqueness, u = u in f2. 

The differentiability properties of u close to M e indicated in Theorems 
1.1.1 and 1.1.2 follow immediately. 

2. The function u does not blow up anywhere else. The proof of the 
upper bound of the lifespan (1.1.1) and of Theorem 1.1.2 is already complete. 
It remains to prove the additionnal statement of Theorem 1.1.1 that u blows 
up only at M £ . 

We recall first the interior asymptotic results of [1] (see Theorem 2.4 and 
Section 2.5). We define S(a, uj,t) to be the solution of 

d T S-^(d a S) 2 = 0, S(<t,w,0) = RfUfaw) +eR ( - 2 \a,uj). 
We then set 

u a = x(j-^l}~T/2 S ( r ~ *' w ' r )> ( d t ~ A )^« ~ (d t u a )(d 2 u a ) = J a , 

where x( s ) 1S zero f° r s < 1/2 and is one for s > 3/2. It turns out that this 
approximate solution satisfies 

(5.2.1) \d^Ju-u a )(.,t)\ < C a e 2 ^\lne\ 
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for r close to to and 

(5.2.2) \d^Ja(.,t)\ < C a e 5 \lne\, \V 2 u a \ < Ce 2 

if the norms are taken only in the "interior" domain r — t < —Co, and t <T e . 

The following lemma explains how this approximate solution can be glued 
together with the exact solution u to yield a new approximate solution u a . 

Lemma V.2. There exists an approximate solution u a with the same 
properties (5.2.1) and (5.2.2) as u a , and moreover u a = u for r — t > —Co- 

Proof, a. Consider G in a strip 

-Ci <r-t< -C 2 

not containing the blowup point. We have 

d a (drG - \{d a G) 2 ) = e 2 [^- + B]{$- 1 ), 

where A and B are smooth. Hence, by integration, we obtain 

d T G-^(d <T G) 2 =e 2 C(^- 1 ), 

with C smooth. On the other hand, the trace on {r = tq} of S and G differs 
by a smooth function which is 0(e 2 lne). Hence G — S = 0{e 2 lne). 
b. If we now set 

Ua = X( r ~ t)u a + (1 - X( r ~ t))u 

for an appropriate %, we can easily check the L 2 estimates of the traces and 
of J a . □ 

We now set u = u a + u and estimate ii using standard energy inequalities 
between tq and f £ . When doing so, we need only control V 2 u a on the support 
of ii, that is, where it is bounded. It follows that V 2 ii, and more generally all 
derivatives of n, are also bounded, which completes the proof. 
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